Determination of neutrino mass hierarchy by 21 cm line and CMB B-mode 

polarization observations 
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We focus on the ongoing and future observations for both the 21 cm line and the CMB B-mode 
polarization produced by a CMB lensing, and study their sensitivities to the effective number of neu- 
trino species, the total neutrino mass, and the neutrino mass hierarchy. We find that combining the 
CMB observations with future square kilometer arrays optimized for 21 cm line such as Omniscope 
can determine the neutrino mass hierarchy at 2<r. We also show that a more feasible combination 
of Planck + Polarbear and SKA can strongly improve errors of the bounds on the total neutrino 
mass and the effective number of neutrino species to be ASm„ ~ 0.12 eV and AN U ~ 0.38 at 2a. 



I. INTRODUCTION 

Since the discoveries of neutrino masses by Super- 
Kamiokande through neutrino oscillation experiments in 
1998, the standard model of particle physics has been 
forced to change to theoretically include the neutrino 
masses. 

So far only mass-squared differences of the neutri- 
nos have been measured by neutrino oscillation exper- 
iments, which are reported to be Amji = rn\ — m\ — 
7-59±~oi? x 10" 5 eV 2 and Am 2 2 = m§ - m\ = 
2.43±o;i3 x lCT 3 eV 2 @|. However, absolute values and 
their hierarchical structure (normal or inverted) have not 
been obtained yet although information for them is in- 
dispensable to build such new particle physics models. 

In particle physics, some new ideas and new future 
experiments based on those ideas have been proposed 
to observe the absolute values and/or the hierarchy of 
neutrino masses, e.g., through tritium beta decay in KA- 
TRIN experiment @ , neutrinoless double-beta decay [U , 
atmospheric neutrinos in the proposed iron calorimeter 
at INO [5j,j6J] and the upgrade of the IceCube detector 
(PINGU) pf, and long-baseline oscillation experiments, 
e.g., NOj/A [|, J-PARC to Korea (T2KK) [HGjJ and Oki 
island (T2KO) [TO], and CERN to Super-Kamiokande 
with high energy (5 GeV) neutrino beam (l2l |. 

On the other hand, such nonzero neutrino masses af- 
fect cosmology significantly because relativistic neutrinos 
prohibit the perturbation from evolving, due to following 
two reasons. First of all, in general relativity, the density 
perturbation of a relativistic particle can hardly evolve 
at all before it becomes nonrelativistic. Second a rela- 
tivistic neutrino erases its own density perturbation up 
to a horizon scale through its free streaming at every cos- 
mic time. By measuring spectra of density perturbations 
by using observations of cosmic microwave background 
(CMB) anisotropics and large-scale structur e (LSS) , we 
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could constrain the total neutrino mass E m v [K 
the effective number of neutrino species N v \2u 
l3lj ]. So far the robust upper bound on E m„ has been 
obtained to be Em„ < 0.62 eV (95 % C.L.) (see Ref. [H 
and references therein) by these cosmological observa- 



tions. For forecasts by future CMB observations, see also 
Refs. H3,[l|. 

In addition, by observing power spectrum of cosmo- 
logical 21 cm radiation fluctuation, we will be able to 
obtain independent useful information for the neutrino 
masses [3J-|37|. That is because the 21 cm radiation 
is emitted (1) long after the CMB epoch (at a redshift 
z <C 10 3 ) and (2) before an onset of the LSS forma- 
tion. The former condition (1) gives us information on 
smaller neutrino mass ( < 0.1 eV). The latter condition 
(2) means we can treat only a linear regime of the matter 
perturbation, which can be analytically calculated unlike 
the LSS case. 

In actual analyses, it should be essential that we com- 
bine data of the 21 cm with that of the CMB observations 
because they complementary constrain cosmological pa- 
rameter spaces each other. Leaving aside minded neu- 
trino parameters, for example, the former is quite sen- 
sitive to the dark energy density, but the latter is rela- 
tively insensitive to it. On the other hand, the former 
has only a mild sensitivity to the normalization of the 
matter perturbation, but the latter has an obvious sensi- 
tivity to it by definition. In pioneering works by [36| , the 
authors tried to constrain the neutrino mass hierarchy 
by combining Planck satellite with future 21 cm obser- 
vations in case of relatively degenerate neutrino masses 
T,m v ~ 0.3 eV. 

Here, however, we additionally include analyses of the 
CMB B-mode polarization produced by a CMB lensing. 
This gives us more detailed information on the matter 
power spectrum at later epochs, which means it has bet- 
ter sensitivities for smaller neutrino masses down to < 
0.1 eV. That is essential to distinguish the normal hi- 
erarchy from the inverted one. In particular we adopt 
ongoing and future CMB observations such as Polar- 
bear and CMBPol, which have much better sensitivities 
to the B-mode. For ongoing and future 21 cm observa- 
tions, we adopt MWA, SKA and Omniscope experiments. 
We forecast possible allowed parameter regions for both 
neutrino masses and effective number of neutrino species 
when we use the above-mentioned ongoing and future 
observations of the 21 cm and the CMB. In particu- 
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lar we propose a nice combination of neutrino masses, 
r„ = (m 3 — mi)/£m„ to make the mass hierarchy bring 
to light as is explained in the text. 



II. 21 CM RADIATION 

Here we briefly review basic methods to use the 21 cm 
line observations as a cosmological probe. For further 
details, we refer readers to Refs. [33.l38|. 



A. Power spectrum of 21 cm radiation 

The 21 cm line of the neutral hydrogen atom is emitted 
by hyperfine splitting of the IS ground state due to an 
interaction of magnetic moments of proton and electron. 
Spin temperature T$ of neutral hydrogen gas is defined 
through a ratio between number densities of hydrogen 
atom in the IS triplet and IS singlet levels, nx/riQ = 
(gi/go)exp(-T*/Ts). where T* = hc/k B X 21 = 0.068 K 
with A2i(— 21 cm) being the wave length of the 21 cm 
line at a rest frame, and gi/go — 3 is the ratio of spin 
degeneracy factors of the two levels. A difference be- 
tween the observed 21 cm line brightness temperature at 
redshift z and the CMB temperature Tcmb is given by 



T 6 (x)« 27x m (l + 6 b ) 



0.15 1 



0.023/ \Q m h 2 10 
( T s -Tcmb \ fH(z)/(l + z) 



dvu/dn 



1/2 



mK, (1) 



where cchi is the neutral fraction of hydrogen, 5b is the 
hydrogen density fluctuation, and dv\\/dr\\ is the gradient 
of the proper velocity along the line of sight due to both 
the Hubble expansion and the peculiar velocity 

In general, T& is sensitive to details of intergalactic 
medium (IGM). However, with a few reasonable assump- 
tions we can omit this dependence [39144 1| . At an epoch 
of reionization (EOR) long after star formation begins, 
X-ray background produced by early stellar remnants has 
heated the IGM. Therefore a gas kinetic temperature Tk 
could be much higher than the CMB temperature Tcmb • 
Furthermore the star formation produces a large back- 
ground of Lya photons sufficient to couple Ts to T# via 
the Wouthuysen-Field effect [IH, Ej|. In this scenario, we 
are justified in taking Tcmb *C Tk ~ Ts at z < 10, so 
that Tb does not depend on Ts- 

In addition, we adopt following assumptions for the 
EOR in the same manner as [36J, |44| . If the IGM is fully 
neutral, fluctuations of the 21 cm radiation arise only 
from density fluctuations. In this limit, we can write the 
power spectrum of the 21 cm line brightness fluctuation 
P 2 i(k)as[H 



Here P 2 i(k,z) is defined by (5T b (k)5T* (k 1 )) = 
(2n) 3 5 3 {k - k')P 2 i(k), where 5T b = T b - T b is the de- 
viation from a spatially averaged brightness temperature 
Tb, Pgs is the matter power spectrum, and /i = k n is the 
cosine of the angle between the wave number k and the 
line of sight. In principle, T, can be calculated although 
it depends on the unknown ionization and thermal his- 
tory. Therefore we treat Tb as a free parameter to be 
measured. 

The power spectrum P2i(k, z) and the comoving wave 
number k are not directly measured by the observations 
of 21 cm radiation [45|,[46[. Instead, here we define u as 
the Fourier dual of © = ^e, + Ojij + A/e^, where 9i and 
9j determine an angular location on the sky plane and 
A/ shows the frequency difference from the central red- 
shift of a z bin. The vector u and its function P2i(u, z) 
are directly measured by the observations. Relationships 
between u = u,ej + Ujij + u\\ik and k are represented by 
= Uiii + ujij = dyt(z)kj_ = 27rL/A, and it|| = y(z)k\\. 
Here "_L" denotes the vector component perpendicular 
to the line of sight. "||" denotes the component in the 
line of sight, ^(z) is the comoving angular diameter 
distance to a given redshift. y(z) = Aai(l + z) 2 /H(z) 
means the conversion factor between comoving distance 
intervals and frequency intervals A/. L is the baseline 
vector of an interferometer. A = A2i(l + z) denotes 
the observed wave length of the redshifted 21 cm line. 
In u space, the power spectrum P 2 i(u, z) is defined by 
(5T b (u)5T b *(u')) = (2tt) 3 5 3 (u- u')P 2 i(u). Then, the re- 
lation between P2i(u, z) and P 2 i(k, z) is given by 



P 2 i(u,z) 



d 2 A {z)y{z) 



P 2 i(k,z). 



(3) 



P 2 i(k,z)=T b 2 (z)(l + M 2 ) P SS (k,z). 



(2) 



We perform our analyses in terms of P2i(u, z) since this 
quantity is directly measurable without any cosmolog- 
ical assumptions. For methods of foreground removals, 
see also recent discussions about independent component 
analysis (ICA) algorithm, FastICA [47| which will be de- 
veloped in terms of the ongoing LOFAR observation [48[ . 



B. Effects of neutrino masses on power spectrum 

The massive neutrinos affect the growth of the matter 
density perturbation mainly due to following two physical 
mechanisms. [49J. First of all, a massive neutrino Vi (even 
with its light mass m Vi < 0.3 eV ) becomes nonrelativistic 
at T ~ m Vi and has contributed to the energy density 
of cold dark matter (CDM), which changes the matter- 
radiation equality epoch and has changed an expansion 
rate of the universe since that time. When we fix the total 
mass of neutrinos T,m u (< 0.3 eV), only the latter effect 
is effective. Second, the matter density perturbations on 
small scales can be suppressed due to the neutrinos' free- 
streaming. As long as neutrinos are relativistic, they 
travel at speed of light, and their free-streaming scales 
are approximately equal to the Hubble horizon. Then 
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the free-streaming effect erases their own perturbations 
within such scales. 

Compared with the standard ACDM models where 
three massless active neutrinos are assumed, we will con- 
sider two more freedoms. First one is an introduction 
of the effective number of neutrino species N v , which 
counts generations of relativistic neutrinos before the 
matter-radiation equality epoch and should not be equal 
to three. Second one is the neutrino mass hierarchy. It 
is clear that a change of N v affects the epoch of matter- 
radiation equality. On the other hand, the neutrino mass 
hierarchy affects both the free-streaming scales and the 
expansion rate as was mentioned above |50l ] . In terms of 
the observations of the 21 cm signal, the minimum cut- 
off of the wave number is giv6n by k m [ n — 2w/(yB) — 
6 x 10 _2 /iMpc~ 1 (see Sec III CP while the wave num- 
ber corresponding to the neutrino free-streaming scale 
is fcfrcc < 10 _2 /iMpc _1 . Therefore the main feature of 
the modification of the matter density fluctuation due to 
the change of the mass hierarchy comes from the mod- 
ification of the cosmic expansion when we fix the total 
matter density at the present time. 



C. Forecasting methods and interferometers 

Here we summarize future observations of the 21 cm 
signals emitted at the EOR. We also provide a brief re- 
view of the Fisher matrix formalism for the 21 cm obser- 
vations. We consider MWA [5l|, SKA [EH and Omnis- 
cope [HI for future observations. The summary of the 
detailed specifications is listed in Table |U Each interfer- 
ometer has its own different noise power spectrum, 



P N (u ± ,z) 



^(z)Tsys(z) 
A e (z) 



1 



t n(u±) ' 



(4) 



which affects sensitivities to the 21 cm signals. Here 
T sys ~ 280[(1 + z)/7.4] 2 3 K is the system tempera- 
ture [H3], to is the total observation time, and A e is the 
effective collecting area of each antenna tile. The effect 
of the configuration of the antennae is encoded in the 
number density of baseline n(u±). In order to calculate 
n(wj_), we have to assume a realization of antenna den- 
sity profiles for each interferometer. For MWA, we take 
500 antennae distributed with a filled nucleus of radius 
20 m surrounded by the remainder of the antennae dis- 
tributed with an r~ 2 antenna density profile out to 750 
m [5|. For SKA, we distribute 20% of a total of 5000 
antennae within a 1 km radius and take the antennae 
distributed with a nucleus surrounded by an r~ 2 antenna 
density profile in the same way as those of MWA. These 
antennae are surrounded by a further 30% of the total 
antennae in a uniform density annulus of outer radius 6 
km [3(| . The remainder of the antennae is distributed at 
larger distances sparsely to be useful for power spectrum 
measurements. Finally, we consider Omniscope that is 
a future square-kilometer collecting area array optimized 



for observations of the 21 cm signal. In case of Omnis- 
cope, we take all of antennae distributed with a filled 
nucleus according to [45| . 

To forecast 1 a errors of cosmological parameters, we 
use the Fisher matrix formalism [56]. For the observa- 
tions of the 21 cm signal, the Fisher matrix for cosmo- 
logical parameters pi is expressed as (57j 



:i21cm 



E 

pixels 



l 



[<SP 2 i(u)] 2 V dp. 



dP 21 (u)\ fdP 21 (u) 



dp 3 



(5) 



where we sum only over half the Fourier space. The 
Fisher matrix determines the errors of the parameter pi 
to be 



A Pl >J{F 



(6) 



The error of the power spectrum measurement SP2i(u) 
in a pixel at u consists of a sum of the sample variance 
and the thermal detector noise. It is expressed as 



dP 2 i(u) = 



Nl /2 



where N c — 2irk±Ak±Ak\\V(z)/ (2ir) 3 is the number of 
independent modes in an annulus summing over the az- 
imuthal angle, V(z) = dA(z) 2 y(z)BxFOV is the survey 
volume, B is the bandwidth, and FOV (ps A 2 /A e ) denotes 
the field of view of the interferometer. For each exper- 
iment, we take account of the presence of foregrounds 
and adopt a cutoff at 2n/(yB) < ku [57?] . We also take 
a maximum value of k to be fc max = 3/iMpc -1 beyond 
which nonlinear effects become important and exclude 
all information for A; max < k. 

For each experiment, we assume a specific redshift 
range as follows 0]. We consider Ly-a forests in ab- 
sorption spectra of quasars and assume that reionization 
occurred sharply at z = 7.5. For an upper limit on the 
accessible redshift range, we take it to be z < 10 be- 
cause of increasing foregrounds and uncertainty in the 
spin temperature at higher redshifts. For the above rea- 
sons, we assume that the observed redshift range of EOR 
is 7.8 — 10.2. Only for MWA, we assume a single redshift 
slice centered at z = 8. 

When we calculate the Fisher matrix, we choose 
the following basic set of cosmological parameters: 
the energy density of matter f2 m /i 2 , baryon flbh 2 , 
dark energy f^A, the scalar spectral index n s , the 
scalar fluctuation amplitude A s (the pivot scale is 
taken to be fc p i vot = 0.002 Mpc~ ), the reionization 
optical depth r, Helium fraction Yne, and the to- 
tal neutrino mass £m„ = mi + mi + m^. Fidu- 
cial values of these parameters (except for Sm„) 
are adopted to be (Q m h 2 , flbh 2 , SIa, n s , A s , t, Yhc) = 
(0.147,0.023,0.7,0.95,24 x lO" 10 , 0.1, 0.24). We set a 
range of the fiducial value of Em„ to be £to„ = 0.05 — 
0.3 eV. Besides these parameters, the brightness tem- 
perature of 21 cm radiation Tb(z) can be taken as a free 
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Experiment N an t 


A e (z = 8 


) ^min 


m ax 


FOV 


z 




[m 2 ] 


[m] 


[km] 


[deg 2 ] 




MWA 500 


14 


4 


1.5 


tt16 2 


7.8-8.2 


SKA 5000 


120 


10 


5 


7r5.6 2 


7.8-10.2 


Omniscope 10 6 


1 


1 


1 


2.1 x 10 4 


7.8-10.2 



TABLE I: Specifications for each interferometers. I/ m i n 
(L max ) is the minimum (maximum) baseline. For MWA, we 
assume a single redshift slice centered at z = 8. For SKA 
and Omniscope, the observed redshift range is z — 7.8 — 10.3, 
and we divide the range into five redshift slices with thickness 
Az ~ 0.5. For each experiment, bandwidth is B — 8 MHz, 
and we assumed observations for 8000 h on two places in the 
sky. We assume that the effective collecting area A e is pro- 
portional to A 2 for MWA and SKA. For Omniscope, both A e 
and FOV are fixed. 



parameter. In this study, we adopt the fiducial values 
of f b (z) to be (f 6 (8),f 6 (8.5),f 6 (9) ) Tj,(9.5),f 6 (10)) = 
(26, 26, 27, 27, 28) in units of mK. 

Additionally, we separately study following two cases: 



Experiment 


V 

[GHz] 


Att 
[ M K-'] 


App 
\isK-'\ 


Cfwhm 
['] 


r 

J sky 


Planck [32J 


70 


137 


195 


14 


0.65 




100 


64.6 


104 


9.5 


0.65 




143 


42.6 


80.9 


7.1 


0.65 


POLARBEAR [64] 


150 




8 


3.5 


0.017 


CMBPol (EPIC-2m) [65] 


70 


2.96 


4.19 


11 


0.65 




100 


2.29 


3.24 


8 


0.65 




150 


2.21 


3.13 


5 


0.65 



TABLE II: Experimental specifications of Planck, Polar- 
bear and CMBPol assumed in this study. Here v is the 
observation frequency, Att is the temperature sensitivity per 
1' x 1' pixel, App is the polarization sensitivity per 1' x 1' pixel, 
#fwhm is the angular resolution defined as the full width at 
half-maximum, and / s k y is the observed fraction of the sky. 
We use ^ max = 2000 for Polarbear, and ^ max = 2500 for 
Planck and CMBPol. 



III. CMB 



(A) Effective number of neutrino species 



CMB and neutrino 



We add one more parameter of the effective number 
of neutrino species N v to the fiducial set of the param- 
eters. The fiducial value of this parameter is set to be 
N v = 3.04. In this analysis, we assumed three species of 
massive neutrinos + an extra relativistic component. 

(B) Neutrino mass hierarchy 

In a cosmological context, many different parameter- 
izations of the mass hierarchy have been proposed [58l — 
[6l| . We adopt r v = (rrc.3 — mi)/£m„ (6l| as an additional 
parameter to nicely discriminate the true neutrino-mass 
hierarchy pattern from the other between the normal and 
inverted hierarchies. The normal and inverted mass hi- 
erarchies mean mi < to 2 <C TO3 and m.3 <C m\ < m-z, 
respectively. We add r v to the fiducial set of the pa- 
rameters. r„ becomes positive (negative) for the nor- 
mal (inverted) hierarchy. It should be a remarkably nice 
point that the difference between tv's of these two hier- 
archies becomes larger as the total mass 'Sm u becomes 
smaller. Therefore r v is particularly useful for distin- 
guishing the mass hierarchy. In Fig [5] we plot behaviors 
of r v as a function of £m„. Note that there is a lowest 
value of Em„ which depends on a type of the hierarchies 
by the neutrino oscillation experiments, i.e., ~0.1 eV for 
the inverted hierarchy and ~0.05 eV for the normal hi- 
erarchy. Therefore, if we could obtain a clear constraint 
like 0.05 eV < £m„ < 0.10 eV, the hierarchy should be 
obviously normal without any ambiguities. As will be 
shown later, however, we can discriminate the hierarchy 
even when 0.10 eV < Y,m v . 



CMB power spectra are sensitive to neutrino masses. 
There are three effects that provide detectable signals for 
the neutrino masses: (1) the transition from relativistic 
neutrino to nonrelativistic one, (2) smoothing of the mat- 
ter perturbation by its free-streaming in small scales, and 
(3) variation of lcnscd CMB power spectra. Future CMB 
experiments are expected to set stringent constraints on 
the sum of neutrino masses [4!| [62[. In particular, the 
last effect is unique in the CMB B-mode polarization pro- 
duced by a CMB lensing. Here we propose to combine the 
CMB experiments with the 21 cm line observations. As 
we will see in Section HVt the combined approach resolves 
degeneracy among some key cosmological parameters and 
is more powerful than individual CMB measurements. In 
addition, it is notable that we are able to detect the ef- 
fective number of neutrino species [63] and determine the 
neutrino mass hierarchy. 



B. Sensitivity and Analysis of the CMB 
experiments 

In this study, we choose Planck [H|, Polarbear [13] 
and CMBPol 65] as examples of CMB experiments. Ex- 
perimental specifications we assumed are summarized in 
Table El 

In our analysis for the CMB, we also take the same 
fiducial model (tt m h 2 , flh 2 , n s , A s , r, Yhc) a s that of 
the 21 cm line experiments (see previous section). We 
evaluate errors of cosmological parameters by using 
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Fisher matrix, which is given by [56 



riCMB \ " (21 + 1) 
*" = > ^ /sky 



x Trace 



dpi 



(8) 



Here is a covariance matrix constructed by using 
CMB power spectra Cf (X = TT, EE, TE), deflection an- 
gle spectrum C dd , cross correlation between temperature 
and deflection angle Cj d , and noise power spectra Nf- 
and Nf d , where C dd is calculated by a lensing poten- 
tial [68[ and is related with Cf B 1 . We compute N dd 
by using a public code FUTURCMB [H which adopts 
the quadratic estimator [68]. In this algorithm, N dd is 
reconstructed by Nj (Y = TT, EE, BB). The covariance 
matrix in the Fisher matrix is expressed as 



G, 



'CJ T +NJ T Cj E 

Cj E C EE + N EE 
Cj d C dd + 



Td 



(9) 



N dd 



where Nj is expressed by using both a beam size 
(v) = 6>fwhm(^)/ v81n2 and an instrumental sen- 
sitivity Ay(^) to be 



n7 = 



E 



i 



Ni (y) 



where 



Nj{v) = A\{u) exp [1(1 + l)aL am (^)] 



(10) 



(11) 



For Aee(i / ) and A^(y), we commonly use App(v) listed 
in Table HU N dd is calculated by iV/ T , N EE , and A^ BB . 

In case of Planck and Polarbear, we combine both 
the experiments, and assume that the 1.7% region of the 
whole sky is observed by both the experiments, and the 
remaining 63.3%(= 65% — 1.7%) region is observed by 
Planck only. Therefore we evaluate a total Fisher matrix 
pCMB ^ summ ing the two Fisher matrices, 



riPlanck 



(/sky = 0.633) 



pPlanck+PB 



(/sky =0.017), 



(12) 



where F planck is the Fisher matrix of the region observed 
by Planck only and F planck+PB is that by both Planck 
and Polarbear. 

In addition, we calculate noise power spectra 
i Pbmck+PB of the CMB polarization ( Y = EE or BB) 

i bmck+PB w j^ n ^he following operation. 



in F t 



1 By performing a public code HALOFIT [66, 67], we have checked 
that modifications by including nonlinear effects for evolutions of 
the matter power spectrum are much smaller than typical errors 
in our analyses and negligible for parameter fittings. 





Sm„ [eV] 


N v 


Fiducial value 


0.05 


3.04 


Planck + Polarbear 0.146 


0.282 


+ MWA 


0.114 


0.240 


+ SKA 


0.0592 


0.189 


+ Omniscope 


0.0226 


0.0753 


CMBPol 


0.0538 


0.0929 


+ SKA 


0.0276 


0.0827 


+ Omniscope 


0.0131 


0.0438 



TABLE III: 1-a e xperimen tal uncertainties of Era„ and N v 
defined by Apt =yj (F^ 1 )^. 



(1) 2 < I < 25, 2000 <£< 2500 



^yY,Planck+PB _ , Planck 



(13) 



Y,Pla„ck +1/iv Y,PB r l (14) 



(2) 25 < £ < 2000 

^Y.Pfanck+PB = ^ 

Since Polarbear observes only CMB polarizations, 

, , , , , atTT, Planck+PB . 

the temperature noise power spectrum l\ e is 

1 i AfTT, Planck 

equal to Nf 

In order to combine the CMB experiments with the 
21 cm line experiments, we calculate the combined fisher 
matrix to be 



p21cm+CMB ^ pCMB , p21cm 



(15) 



Here we did not use information for a possible correlation 
between fluctuations of the 21 cm and the CMB. 



IV. RESULTS 

In this section, we numerically evaluate how we can de- 
termine (A) the effective number of neutrino species, and 
(B) the neutrino mass hierarchy, by combining the 21 cm 
line observations (MWA, SKA, or Omniscope) with the 
CMB experiments (Planck + Polarbear, or CMBPol). 
To obtain Fisher matrices we use CAMB [6^, [67} for cal- 
culations of CMB anisotropics C/ and matter power spec- 
tra Pss- 



A. Constraints on N v 

In FigfTJ we plot contours of 90%confidencc levels 
(C.L.) forecasts in Em^-TV^ plane. The fiducial values 
of neutrino parameters, N v and Sto„, are taken to be 
N v = 3.04 and T,m u = 0.1 eV (upper two panels), which 
corresponds to the lowest value of the inverted hierar- 
chy model, or £m„ = 0.05 eV (lower two panels), which 
corresponds to the lowest value of the normal hierarchy 
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2.5 





Planck+POLARBEAR 




+ MWA 




+ SKA 




" + Omniscope 

\ \ ^ 

\ 1 x 




/ / ' 
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0.1 0.2 0.3 0.4 0.5 0.6 
Em v [eV] 




Planck+POLARBEAR - 

+ MWA- 
+ SKA- 
+ Omniscope - 



0.1 0.2 0.3 0.4 0.5 0.6 
2m v [eV] 



3.4 



3.2 



2.8 



2.6 



2.6 



CMBPol - 

+ SKA - 
■ Omniscope - 




0.05 0.1 0.15 0.2 
Sm v [eV] 




CMBPol ■ 

+ SKA ■ 
+ Omniscope - 



0.05 0.1 0.15 
T.m v [eV] 



0.2 



FIG. 1: Contours of 90% C.L. forecasts in Y>m v -N v plane, by adopting Planck + Polarbear + each 21 cm experiment (left 
two panels), or CMBPol + each 21 cm experiment (right two panels). Fiducial values of neutrino parameters, N v and Em„, are 
taken to be N v — 3.04 and Em„ =0.1 eV (for upper two panels) or Sm„ = 0.05 eV (for lower two panels). The dashed line 
means the constraint obtained by only a CMB observation such as Planck + Polarbear alone (left two panels), or CMBPol 
alone (right two panels). The severer constrains are obtained by combining the CMB with a 21 cm observation such as MWA 
(outer solid, only for left panels), SKA (middle solid), and Omniscope (inner solid), respectively. 



model. Adding the 21 cm experiments to the CMB ex- 
periment, we see that there is a substantial improvement 
for the sensitivities to Yim v and N v . That is because 
several parameter degeneracies are broken by those com- 
binations, e.g., in particular Xj, and As were completely 
degenerate only in 21 cm line measurements. Therefore 
it is essential to add the CMB to the 21 cm experiment 
to be vital for breaking those parameter degeneracies. 

If each CMB experiment is combined with SKA or 
Omniscope, the corresponding constraint can be signifi- 
cantly improved. We showed numerical values of those 
errors in Table IIIII in case that the fiducial values are 
taken to be N v = 3.04 and Y,ra u = 0.05 eV. On the 
other hand, comparing those values with the current best 
bounds for £m„ + N v model, which give £m„ < 0.89 
eV and N v = 4.47t^ obtained by CMB (WMAP) + 
HST (Hubble Space Telescope) + BAO [Hj], we find that 
the ongoing and future 21 cm line -I- the CMB obser- 
vation will be able to constrain those parameters much 



more severely. 

The case of £m„ = 0.1 eV to be fiducial (upper two 
panels) corresponds to the lowest value for the inverted 
hierarchy when we use oscillation data. Then it is notable 
that CMBPol + SKA can detect the nonzero neutrino 
mass. Of course, Planck + Polarbear + Omniscope 
and CMBPol + Omniscope can obviously do the same 
job. 

On the other hand, the case of Sm„ = 0.05 eV to 
be fiducial (lower two panels), which corresponds to the 
lowest value for the normal hierarchy, only Planck + Po- 
larbear + Omniscope or CMBPol + Omniscope can 
detect the nonzero neutrino mass. 



B. Constraints on neutrino mass hierarchy 

Next we discuss if we will be able to determine the 
neutrino mass hierarchies by using those ongoing and fu- 
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FIG. 2: Forecasts of 2a errors on r v = [mz — mi)/Sm„ constrained by both the 21 cm and the CMB observations in case 
of the inverted hierarchy to be fiducial (left), and the normal hierarchy to be fiducial (right). The constrains are obtained 
by combining Omniscope with Planck + Polarbear (thick dashed lines), and Omniscope with CMBPol (thick solid lines), 
respectively. Allowed parameters on r„ by neutrino oscillation experiments are plotted as two bands for the inverted and the 
normal hierarchies, respectively (the name of each hierarchy is written in the close vicinity of the line). The solid line inside 
the band is the fiducial value of r„ as a function of Em„. 



ture 21 cm and CMB observations. In Fig. [2] we plot 
2a errors of the parameter r„ = (7713 — mi)/Em„ con- 
strained by both the 21 cm and the CMB observations 
in case of the inverted hierarchy to be fiducial (left), and 
the normal hierarchy to be fiducial (right). It is notable 
that the difference between r„ 's of these two hierarchies 
becomes larger as the total mass £m„ becomes smaller. 
Therefore, r„ is quite useful to distinguish a true mass 
hierarchy from the other. Allowed parameters on r v by 
neutrino oscillation experiments are plotted as two bands 
for the inverted and the normal hierarchies, respectively. 
The thin solid lines inside the bands are the experimen- 
tal mean values by oscillations, one of which is taken 
to be a corresponding fiducial value of r v as a function 
of £m„ in each analysis. The constrains are obtained by 
combining Omniscope with Planck + Polarbear (thick 
dashed lines) and Omniscope with CMBPol (thick solid 
lines), respectively. 

For 0.3 eV < YjTJijy the mass eigenvalues ttl\, vri2, and 
m3 are almost degenerate. Therefore the difference be- 
tween two hierarchies has little influence on the matter 
power spectrum. Therefore the constraints on r v are 
significantly weak compared with the difference between 
them, and then we cannot distinguish the true hierarchy 
from the other. 

On the other hand however, the difference increases as 
£to„ decreases down to m v ~ 0.1 eV. By using this prop- 
erty, the CMB (Planck + Polarbear or CMBPol) + the 
21 cm (Omniscope) observations can constrain the neu- 
trino mass hierarchy severely. Typically errors of Era^ at 
around Y>m v = 0.1 eV are given by A£m„ = 0.0087 eV 
for Planck + Polarbear, and AS™,, = 0.0069 eV for 
CMBPol at la, respectively. Therefore, the error of the 
x-axis is negligible compared with that of y-axis. In 



Fig. |21 we plot contours of 90% C.L. in T,m v —r v plane in 
order to show errors of Y>m v along the cc-axis for typical 
fiducial values. As is clearly shown in Fig. [2] actually 
those combinations of the observations will be able to 
determine the neutrino mass hierarchy to be inverted or 
normal for Em,, < 0.13 eV or Em„ < 0.1 eV at 90 % C.L., 
respectively. Although the determination is possible only 
at around Em„ < O(0. 1) eV, those results should be rea- 
sonable. That is because a precise discrimination of the 
mass hierarchy itself may have no meaning if the masses 
are highly degenerate, i.e., if £m„ 3> 0.1 — 0.3 eV. 

Once a clear signature Em,, <C 0.1 eV were determined 
by observations or experiments, it should be obvious that 
the hierarchy must be normal without any ambiguities. 
On the other hand if the hierarchy were inverted, we 
could not determine it only by using Em„. However, it 
is remarkable that our method is quite useful because we 
can discriminate the hierarchy from the other even if the 
fiducial values were £m„ > 0.1 eV for both the normal 
and inverted cases. This is clearly shown in Fig [3] In 
case that a fiducial value of Y>m u is taken to be the low- 
est values in neutrino oscillation experiments, the upper 
left (right) figure indicates that even CMBP0I+ SKA can 
discriminate the inverted (normal) mass hierarchy from 
the normal (inverted) one. 

V. CONCLUSIONS 

We have studied how we can constrain effective num- 
ber of neutrino species N v , total neutrino masses Em„, 
and neutrino mass hierarchy by using the 21 cm obser- 
vations (MWA, SKA, and Omniscope) and the CMB ob- 
servations (Planck, Polarbear, and CMBPol). It is 
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FIG. 3: Contours of 90% C.L. forecasts in £m„-r„ plane, by adopting CMBPol + each 21 cm experiment. Fiducial value 
of Em„ and the mass hierarchy (diagonal cross) are taken to be: Em„ =0.1 eV and the inverted (for left upper panel), 
Em„ = 0.12 eV and the inverted (for left lower panel) , Sm„ = 0.06 eV and the normal (for right upper panel), Em„ =0.1 eV 
and the normal (for right lower panel). The short dashed lines mean the constraints obtained by only a CMBPol observation, 
and the long dashed line means the one by a SKA + CMBPol observation for 24000 h on four places in the sky. The outer 
(inner) solid line means combining the CMBPol with SKA (Omniscope). Allowed parameters on r v by neutrino oscillation 
experiments are plotted in the same way of Fig. [2] 



essential to combine the 21 cm with the CMB B-mode 
polarization produced by a CMB lensing to break var- 
ious degeneracies in cosmological parameters when we 
perform multiple-parameter fittings. 

About the constraints on Yim v -N v plane, for a fidu- 
cial value Y>m u =0.1 eV which corresponds to the low- 
est value in the inverted hierarchy, we have found that 
CMBPol + SKA, Planck + Polarbear + Omniscope 
and CMBPol + Omniscope can detect the nonzero neu- 
trino mass. For a fiducial value T,m u = 0.05 eV, which 
corresponds to the lowest value in the normal hierarchy, 
Planck + Polarbear + Omniscope or CMBPol + Om- 
niscope can detect the nonzero neutrino mass. 

As for the determination of the neutrino mass hier- 
archy, we have proposed a new parameter r v = (7713 — 
mi)/Em„ and studied how to discriminate a true hier- 



archy from the other by constraining r v . As was clearly 
shown in Fig.[5J the combinations of the CMB (Planck + 
Polarbear or CMBPol) + the 21 cm (Omniscope) will 
be able to determine the hierarchy to be inverted or nor- 
mal for £m„ < 0.13 cV or < 0.1 eV at 2a, respectively. 
Furthermore, if the fiducial value of Y>m v is taken to be 
the lowest value in the neutrino oscillation experiments, 
even CMBPol + SKA can determine the mass hierarchy. 

In this study we have taken the simplified model of 
rcionization. In case of more likely detailed modeling of 
reionization [45| . it was pointed out that the constraints 
on cosmological parameters may moderately change at 
~ 10—50 %. Fortunately, this effect is comparatively 
small and should not be fatal to constrain the neutrino 
mass hierarchy in the current analyses. 
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